Vertex Models and Random Labyrinths: Phase Diagrams for Ice- type Vertex Models. 
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We propose a simple geometric recipe for constructing phase diagrams for a general class of vertex models 
obeying the ice rule. The disordered phase maps onto the intersecting loop model which is interesting in its own 
right and is related to several other statistical mechanical models. This mapping is also useful in understanding 
some ordered phases of these vertex models as they correspond to the polymer loop models with cross-links in 
their vulcanised phase. 



I. INTRODUCTION 

Ice-type models were originally introduced in order to de- 
scribe the properties of ice[lj |2J], and have been later gener- 
alised to represent other types of hydrogen-bonded crystals 
01- In ice, the oxygen atoms form a lattice with the coordina- 
tion number of four. Each bond of the lattice contains a single 
hydrogen atom that is shifted from the middle point toward 
one of the neighbouring oxygen atoms. It has been conjec- 
tured, on the basis of local electro-neutrality, that precisely 
two hydrogen atoms are located near each oxygen atom with 
the other two being shifted away from it. This rule is known as 
the ice rule, and it can be graphically represented by placing 
polarisation arrows along the hydrogen bonds: For each site 
on the lattice there are exactly two incoming and two outgoing 
arrows. There are six such possible arrow configurations at 
each site which leads to another common name for this model 
(or rather its two-dimensional square lattice version): the six- 
vertex model. Interestingly, the states of ice-type models can 
be characterised by topological winding numbers and the ex- 
citations are topological defects carrying fractional charge 01, 
hence the recent revival of interest to this type of models in 
the context of topological order and fractionalisation in quan- 
tum systems |5l la 01- The quantum phase transition in the 
quantum square lattice six-vertex model have been argued to 
mimic the physics of the HTSC materials with the li-density 
wave order while their counterparts on the triangular lattice 
twenty-vertex model could be relevant to the physics of 2D 
Wigner crystals 1 8 ] . 

The six-vertex model on the square lattice was solved ex- 
actly by Lieb j3 E El El and Sutherland Q using 
the Bethe Ansatz and also (for some cases) by mapping it 
to a soluble free fermion problem I14ll . For a review see 
OB. Much less is known about the ice-rule models in 
other contexts. For example, in the "next" 2D model, namely 
the twenty-vertex model on the triangular lattice, the Bethe 
Ansatz only works in a few instances when there are special 
relations between the parameters Odd. More impor- 
tantly, there are very few broad-based techniques of general 
applicability with which these sorts of models can be studied. 

In this paper we develop a new approach to establishing 
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and understanding the phase diagrams of vertex models with 
the ice rule (i.e. equal numbers of incoming and outgoing 
arrows at each vertex). We introduce a new class of poly- 
mer models which is closely related to loop models, recently 
a topic of intensive study (see 12011 and references therein). 
While a lot of progress has been achieved in studying loops 
models by the means of exact solutions [21] and conformal 
field theory [22], in the present analysis we will not recurse to 
those methods. The proposed general model is closely related 
to so-called Lorentz lattice gas ll23l l24l EM Erjll or random 
labyrinths 82711 as well as the ice-rule vertex models as will 
be discussed below. Such unified graphical representation on 
which our approach is based, is not entirely new, it has been 
discovered in the context of a loop algorithm developed to 
simulate vertex models 1 281 12^1 . Notwithstanding the compu- 
tational benefits inherent in this representation, here we will 
concentrate on the physical insight into the problem and the 
theoretical benefits which it provides. Indeed, our perspective 
allows for a very intuitive approach to a general class of vertex 
models and opens new possibilities for their rigorous analysis. 
This approach is ideologically similar to that of Fortuin and 
Kasteleyn 1 30] for the case of the Potts model. 

Since a lot is known about the square lattice six-vertex 
model, this is a good reference point from which to start. We 
begin by reviewing a connection between the six-vertex and 
an intersecting loop model and show how the phase diagram 
can be immediately inferred. We then apply our method to a 
particularly interesting case of the twenty-vertex model iden- 
tifying its phase boundaries. A more complete treatment of 
this problem with necessary proofs and additional examples 
shall be presented elsewhere lBlll . 



II. THE SWITCH MODEL 

Let us first introduce a statistical model of objects that we 
call switches. Consider a finite lattice L with an even coor- 
dination number. For simplicity let us discuss the case of a 
homogeneous lattice and, in particular, postpone the consider- 
ation of boundary sites. A switch is a variable associated with 
a lattice site. It is defined as a sorting of all incident edges 
into associated pairs. A more formal definition is as follows: 
Consider a complete graph G(i) whose vertices are the bonds 
of L connected to the lattice site i, A switch otjt (i) is a 
perfect matching in G(z'). Clearly, if 2m is the coordination 
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number, there are (2m — 1)!! possible switches for every lat- 
tice site. On the square lattice, the three possible switches, the 
oc-switch, the p-switch and the y-switch are shown in Fig.^ In 
Ref. 12711 . in the context of random labyrinths, these have been 
referred to as NW and NE mirrors and tunnels respectively.. 
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FIG. 1: The three possible switches for the two-dimensional square 
lattice. We do not distinguish between under- and over-crossings for 
the gamma-switches (we are not concerned with topology for the 
purposes of this paper). 



In general, we shall call the switches the ai -switch, the 
a2-switch etc. We will denote by S a configuration of such 
switches on L. The weight of S , which we denote by Q(s) 
is given as follows: First, there are a priori weights (or activ- 
ities) associated with every type of switches. These will be 
conveniently denoted by 0Ci, OC2 etc. Next, we observe that S 
divides the bonds of the lattice into loops. Indeed, each bond 
is a member of two switches ttjc(i) and OCj/Q') S S which 
pair it with two other bonds etc. We can continue growing 
this path in both directions until the loop is closed, i.e. the 
two "end" bonds are paired by a switch. A loop is thus a cy- 
cle of bonds such that every two adjacent bonds are paired 
up by a corresponding switch in S . We are not distinguish- 
ing the cycles different only by a reversal of the overall order. 
Notice that in general such loop can visit the same site more 
than once and in particular, due to tunnels, can have multi- 
ple self-intersections. We let t(s) denote the number of such 
loops and Ai (5), A2 (s), ■ . .the number of oci-switches, 0C2- 
switches etc. Then Q is given by 
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where p = (2m — 1 ) ! ! is the total number of distinct switches 
and n is a positive real number. If n is an integer, it can be 
thought of as a number of different colours available for every 
loop. The partition function is simply 



z = z ai ...., av =z = £e(5). 



(2) 



that all such fused loops count as one for the purposes of eval- 
uating the corresponding weights in Eq. ([Q. 

As has already been mentioned, the switch model and 
its extension are closely related to several known statistical- 
mechanical models. 

The n — > limit of the switch model describes the Eulerian 
walks or cycles (depending on the boundary conditions), i.e. 
walks that traverse every bond exactly once. Self-intersections 
are allowed via tunnels. If one were to adopt such a model 
for describing a polymer, the weight for the tunnels should 
translate into the rigidity of a polymer molecule. 

The n = 1 case corresponds to Lorentz lattice sas l23l l24l 
I25H26I1 or random labyrinths l27ll . The model in Il27ll on a hy- 
percubic lattice 1r is nearly identical to our extended model, 
with the following difference: instead of cross-links, there are 
normal sites at which a path that is otherwise consistent with 
the switches can randomly, with equal probability, go in any of 
the 2d directions. This is equivalent to constructing a random 
walk on a graph whose vertices are the cross-links and edges 
are the switch-mediated paths between them. The model cor- 
responds to n = 1 since the switches (mirrors and tunnels) as 
well as cross-links (normal sites) are distributed with indepen- 
dent probabilities. 

For arbitrary n, if all switch weights are identical, this is 
just a fully-packed limit of the O(n) loop model l32ll . 

A special integrable case of a model ([TJ with n < 2 corre- 
sponds to a supersymmetric spin chain 1 33 ] 

Finally, the n = 2 case corresponds to the ice-rule vertex 
models, this is the subject of the following section. 



III. RELATION TO THE ICE-TYPE MODELS 

A mapping between n = 2 switch models and and ice-type 
models is easiest seen in the case of the square lattice six- 
vertex model whose possible vertices are shown in Figure |2] 
Each of those vertices is assigned a corresponding energy e, 

*h Hh 4* 

1 2 3 4 5 6 

FIG. 2: The six vertices allowed by the ice rule. 



We also introduce a simple extension of the above switch 
model in which we allow two (or more - depending on the 
lattice) paths at a given site to be "fused" together. We shall 
refer to such objects as cross-links. In general, the cross-links 
may be partial in the sense that not all pairs of bonds adja- 
cent to a site are fused. For a four-coordinated lattice such as 
the 2D square lattice, there naturally can be only a single type 
tying up all four adjacent bonds. Aside from associating the 
additional a priori weight (j) with such cross-links, their intro- 
duction results in the following. If n is an integer, loops that 
have been fused together are now required to have the same 
colour. For an arbitrary n this translates into the requirement 



(i = 1,...,6). In the absence of external electric field, the 
model should be invariant under the reversal of all arrows thus 
requiring £1 = £2, £3 = £4 and £5 = Eg. For the remainder of 
this paper we shall concern ourselves not with the energies as 
such, but rather with the (unnormalised) Boltzmann weights 
Wj exp(— fte,). Conforming to to the usual notations (cf. I15l 
UaO . we shall use a = W\ = w%, b = W3 = W4 and c = W5 = wg. 

The mapping of the switch model to the six-vertex model 
is done by assigning directions to all loops (since there are 
two choices for every loop, the factor of 2 f in Eq. Q can be 
obtained by summing over all such choices). Tracing out the 
switch variables as illustrated in Fig. [2 a) leaves us with the 
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six-vertex configurations [47]. The weights are given by 



a = P + y, b = a + y, c = a + (3 



(3) 
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FIG. 3: (a) Tracing out switches by fusing directed loop segments 
into the ice-type vertices (only types 1, 3 and 5 are shown, types 2, 
4 and 6 are obtained by merely reversing all arrows), (b) The phase 
diagram of the six-vertex model. 

In general, we claim that the partition function for the 
switch model with weights Q where n = 2 is exactly the parti- 
tion function for an ice-rule vertex model with vertex weights 
determined by 



k 
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where A* = 1 if the vertex V, and the switch Sk are consis- 
tent and is otherwise. V, and Sk are consistent if there is 
one incoming and one outgoing arrow for each pairing of the 
switch. (Notice that this automatically preserves the ice rule.) 
The most direct way to see this is to consider a joint measure 
on both a switch configuration S and an arrow configuration 
A with the weight 
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where A(s,Sl) = 1 if S and A are consistent in the above sense 
and otherwise. Now, summing this weight over all possi- 
ble arrow configurations A leads to the the weight of S given 
by Eq. Q while tracing out switches leaves us with ice-type 
arrow configurations where the weight for a single vertex is 
given by Eq. (|4j- 

We now take a closer look at Eqs. (0 and notice that if 
a, p,y > (and the switch model only makes sense for non- 
negative weights), the allowed values of a, b and c lie pre- 
cisely within the entire disordered (DO) region of the six- 
vertex model (Fig. |3jb)). In fact, we argue, this is not a co- 
incidence but a general case: a switch model with all switch 
weights positive always maps onto a DO phase of a vertex 
model 1 48]. But before we present our argument, we should 
discuss the boundary conditions which are defined as follows: 
let us regard our lattice L as actually sitting inside the larger 
(possibly infinite) A. Let 3L denote those bonds that connect 
a site in L with a site outside. A boundary condition is defined 
as a constraint on these bonds' colours (for n integer) or their 



external connectedness (general n). From the perspective of 
the partition function, we note that the paths in L either form 
closed loops or end at the boundary; each boundary bond now 
serves as a "source/sink" for a path. 

There are several rules that can be now invoked: The sim- 
plest two - which are the extreme cases - are (i) all such paths 
count as additional loops in Eq. for the weight of a given 
switch configuration, or (ii) the only loops counted are the in- 
ternal loops. There are also many intermediate sets of rules as 
well as superpositions of various boundary conditions, each 
with their own set of rules. Both free and periodic boundary 
conditions are examples of such superpositions. It is noted 
that none of the above amounts to any internal constraints - 
they are all "external constraints", i.e. all internal switch con- 
figurations S are legitimate. This is in sharp contrast to bound- 
ary conditions that are imposed by placement of arrows. In- 
deed, consider the six-vertex model on a cylinder of girth L. If 
m upward arrows and L — m downward arrows are placed on 
the bottom boundary, then the only allowed switch configura- 
tions are those in which there are \L — 2m\ paths that start at 
the bottom and end at the top. Under such circumstances, it is 
easy to see that boundary conditions are capable of affecting 
the bulk free energy. 

What about the order parameter for the switch model? For 
an integral n, the natural choice is the expectation value for a 
bond at the origin to be of a particular colour. Similarly to the 
Fortuin-Kasteleyn (FK) random cluster model, this translates 
into the probability that a given bond at the origin is connected 
to the boundary; with this definition n needs not remain an in- 
teger. Two bonds are considered connected if and only if there 
exist a switch-mediated path of bonds between them (using 
cross-links is also allowed in constructing such path). In the 
absence of cross-links, this means that the two bonds belong 
to the same loop. If cross-links are allowed, the two bonds 
must belong to the same cluster of fused together loops. 

However, there is an important difference between the or- 
der parameter in the FK model and that in our model. While 
in the integer q FK model, the probability that the origin is 
connected to the boundary is (at least for the wired bound- 
ary conditions) equal to the spontaneous magnetisation of the 
corresponding Potts model II34II . in our model such relation 
between the graphical and the physical order parameters be- 
comes more complicated. To illustrate this point, consider 
again the six-vertex model. Unlike the Ising magnet, it actu- 
ally has two ordered phases, ferroelectric (FE) and antiferro- 
electric (AFE) that are not related to each other by a simple 
gauge transformation. While for each given instance of the 
six-vertex model only one of these ordered phases can be ac- 
cessed by the means of lowering temperature (which we have 
absorbed in the vertex weights), there are two possible order 
parameters with no obvious relation between them. One of 
them, namely the uniform polarisation is discussed in Ap- 
pendix |X] For now it suffices to note that connectedness to 
the boundary is a necessary, if not sufficient condition for any 
type of physical order; the internal loops can be reversed at 
will and thus do not contribute to a long-range order. Along 
the same lines, we note that any two locally defined quantities 
at two locations can be correlated only if the corresponding 
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bonds belong to the same loop. The nature of such correlation 
is trivial in terms of colours but may be much more compli- 
cated in terms of arrows. In this sense, the probability of two 
bonds to belong to the same loop is an upper bound for any 
(truncated) correlation function. 

Let us now present a physical (albeit non-rigorous) argu- 
ment supporting our earlier statement that a n = 2 switch 
model with all switch weights positive is always in a DO 
phase. All we have to show is that less than a finite frac- 
tion of bonds in the interior are connected to the boundary. 
We start by considering a n = 1 model in an upper half-space. 
Because n = 1, we do not need to count loops and may con- 
struct the partition function by "growing" our system upward 
from the boundary layer-by-layer, the choice for each switch 
being completely independent of others. In such a setting it is 
easy to see that the assumption of a finite limit of density of 
(vertical) bonds connected to the bottom is inconsistent with 
the finite probability of terminating any two such strands with 
finite separation between them (by turning paths toward each 
other with the analogs of a- and p-switches and then connect- 
ing them via a finite chain of tunnels) at every given level. 
This argument does not preclude a power-law or even loga- 
rithmic decay of such density away from the boundary, and in 
fact this is precisely what has been observed numerically 12411 . 
What about n = 2, the case of our primary interest? Lacking a 
proo f (the previous part of the argument can be made rigorous 
13111 '). one can nevertheless argue that the n > 1 case always 
favours shorter loops thus strengthening the above conclusion 
ll49ll . Notice that there is nothing particularly 2D about this 
argument, we expect our criterion to work in any dimension- 
ality. 

Let us now turn to the ordered phases. As follows from the 
above discussion, a long-range order can not appear unless 
one of the switch weights vanishes. To see that this condi- 
tion indeed correctly describes the boundary, it is instructive 
to look at the six-vertex model again. Let us start with the FE 
phase boundary, e.g. the one described by b = a + c or, in our 
language, p = 0. Notice, that a line going up can now either 
continue going up passing through y-switches or turn to the 
left via an oc-switch, upon which it will continue to the left 
until the next oc-switch will send it up again. It can never turn 
back on itself! Thus, the total number of loops I in Eq. Q in 
an L x L sample is at most 2L and scales sublinearly with the 
size of the system N = L 2 . As a result, in the thermodynamic 
limit, the factor of 2 in Eq. Q is inconsequential and can be 
dropped. Then the partition function is trivially calculated as 

Z= £ a*f- A = (a+y) N = b N (6) 

se{a,y} N 

meaning that the system becomes fully frozen: ether all ver- 
tices are of type 3 or of type 4. Notice that once again, this 
argument is not specifically 2D and could be easily used in 
higher dimensions, e.g. it generalises the proof by Nagle ll35ll 
of a first order transition in a KDP model on d-dimensional 
tetrahedral lattices. Inside this phase, i.e. when b > a + c, 
we need to extend the model by allowing cross-links so that 
b = oc+y+<|>, a = y, c = a (this is in fact very similar to 
"freezing" a vertex in the loop algorithm of 12811 '). Notice that 



in general a cross-linked model is not automatically equiv- 
alent to an ice-type vertex model. One has to verify that a 
consistent arrow assignment can be done for all cross-linked 
clusters with the cross-links corresponding to the excessive 
vertices. For the FE phase, cross-linking the switch model is 
automatic, but not particularly enlightening since the model 
fully orders already at the transition point. The AFE phase 
presents a more interesting case. At the boundary of the AFE 
phase y = 0, i.e. loops turn at every step and cannot self- 
intersect. Assigning directions renders each such loop consis- 
tent with a perfect staggered AFE order, but since each loop 
is free to choose one of the two possible directions, no over- 
all order results. Moving into the phase, we once again resort 
to cross-links to compensate for excessive c: c = a + P + (|), 
a = P, b = a. Cross-links do not change the geometry of loops 
but rather "vulcanise" them into clusters. It is easy to verify 
that such clusters are always consistent with one of the two 
possible staggered polarisations. What is essential here is the 
fact that in the DO phase the loops are critical, as expected in 
this type of models on general grounds. Therefore such vul- 
canisation leads to percolation of a particular colour or, alter- 
natively, staggered polarisation throughout the sample. What 
is a bit puzzling, is that the AFE transition is actually of the 
infinite order, not of the second order as our naive argument 
would generally suggest; this must be due to the criticality of 
the phase on the DO side but we do not have a clearer picture 
for that in terms of loops. 



IV. TWENTY- VERTEX MODEL ON THE TRIANGULAR 
LATTICE 

So far we dealt with the lattice of coordination number four. 
This is somewhat special since the number of possible vertices 
(up to the reversal of all arrows) allowed by the ice rule is the 
same as the number of possible switches - three. This is not 
generally the case. Let us turn our attention to the triangular 
lattice. Its coordination number is six which means there are 
(3) = 20 distinct ways of arranging tree inward and tree out- 
ward pointing arrows. By requiring that the weights for the 
vertices are arrow-reversal invariant (zero electric field con- 
dition), the number of parameters is reduced to 10. On the 
other hand, there are 5!! = 15 distinct switches. Therefore 
Eqs. (0J form an underdetermined system. So, at a first glance 
the situation appears hopeless: we cannot even uniquely solve 
this system, so how could we expect to deduce the phase 
boundaries? Nevertheless, the proposed criterion stands: as 
long as it is possible to satisfy Eqs. (|4} with positive switch 
weights, the system is in its (critical) DO phase. To illus- 

FIG. 4: The distinct (up to rotations and reflections) vertices and 
switches on the triangular lattice with their weights. The third vertex 
is chiral, with only one chirality shown; both have the same weight. 
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trate this, consider the symmetric version of the twenty-vertex 
model (dubbed the F-model in 11711 ) by requiring that all ver- 
tex weights are invariant under the symmetry group D(,. This 
leaves us with only three distinct types of vertices and five 
types of switches depicted in Fig.0] Eqs. (0} read: 

a = p + 2y+28 + £, 

£ = 2a + 3p + e, (7) 
c = a + p + 3y+S, 

or, equivalently 

a + 2c = Z? + 8y+48, 
/? + 2c = a + 4(a + p + y), (8) 
£ + 3a = 2c + 4(p + 8 + £). 

From the latter form, it is clear that there are no solutions with 
all switch weights positive unless a + 2c > b, b + 2c > a, b + 
3a > 2c; according to our criterion this is the DO region of the 
model. The suspected phase boundaries correspond to a = 
p = y = 0, y = 8=0 and p = 8 = e = 0. These conditions 
have a simple geometric interpretation in terms of loops: in 
the first case, connections of bonds at 60° are excluded, in the 
second, the are no connections at 120° and in the third, there 
are no direct tunnels. We now turn to the suspected ordered 




FIG. 5: (a) The proposed phase diagram for the F-model on the trian- 
gular lattice shown. Along the dashed line the Bethe Ansatz solution 
is available with the known transition point circled, (b) Complete and 
partial cross-links important for understanding the ordered phases II 
and I (see text for details), (c) A zero mode in a state composed 
entirely of c-vertices; the arrows along the selected pentagon can be 
freely reversed. 

phases. The simplest one is the phase II (b > a + 2c); it is 
easy to verify that the introduction of complete cross-links ((j)* 
in Fig. |5j?) corresponding to excessive ^-vertices is consistent 
with the onset of the spontaneous sextipole electric moment 
- a triangular flux phase. This is the only ordered phase of 
the model accessed by the means of exact solutions | 3 Inil 
available along the dashed line in Fig. \5]p, a transition to the 
ordered phase was found to be of the infinite order. 

In the phase I (a > b + 2c), on the other hand, complete 
cross-links are inconsistent with the dipolar order described 
by the proliferation of a-vertices; such cross-link correspond- 
ing to any particular a-vertex cannot be broken into all three 
different 8-switches. Physically, this inconsistency arises 
from the fact that the symmetry of the ordered phase (mod 



complete arrow reversal) is lower than the symmetry of the 
cross-link and switch weights. One possible remedy is to re- 
place complete cross-links with the partial ones (0 X in Fig.[5p) 
breaking such symmetry by hand. The particular choice of 
two out of three possible partial switches prefers a horizontal 
dipolar order: a_ = 28 + £ + 2§ x , a\ = a/ = 28 + £ + (|)x 
where the subscripts for a-vertices indicate directions of their 
middle out-arrows. Now once again, vulcanisation of loops 
leads to the onset of the dipolar order in the vertex model, but 
this is not the original F-model, we have reduced the symme- 
try of the a-vertex weights. We nevertheless expect a simi- 
lar transition to occur in the F-model since the two models 
are identical at the transition line of the modified model and 
their ordered phases have the same broken symmetry 
plausible counter-argument would point at the criticality of 
the DO phase making the above reasoning suspect: moving 
in the "wrong" direction on the phase diagram can easily or- 
der a critical phase while the system might remain disordered 
if we carefully tread in the "right" direction. We argue, this 
is not the case: the disordered phase here actually expected 
to have finite susceptibility to the uniform electric field as is 
the case for the six vertex model 1 15, 36] |51]. Therefore the 
non-analytic onset of the dipolar order in any direction of the 
phase diagram signifies a true phase transition. 

Finally, there is Region III (2c > b + 3a) described by the 
abundance of c-vertices. Cross-linking the loops does not ap- 
pear to work here since unlike in the previous two cases, some 
ground states of system have local zero modes (Fig.|5j;). What 
is even more unsettling, is that in terms of loops a particular 
zero mode depicted in Fig.|5j; contains a direct tunnel, some- 
thing that disappears on the boundary of Region III but seem- 
ingly reappears well inside. We therefore can only speculate 
on the nature (and the existence) of the transition into this 
phase, although the fact that a certain class of loops that could 
freely fluctuate in the DO phase disappears at the boundary 
remains very suggestive. One possibility is the onset of an 
entropically stabilised phase, i.e. an order-by-disorder transi- 
tion; further investigation is needed to clarify this issue. 



V. TWENTY-VERTEX MODEL ON THE CUBIC LATTICE 

Interestingly, there is a one-to-one correspondence between 
the vertices (and corresponding switches) in twenty-vertex 
models on the 2D triangular and 3D cubic lattice (it becomes 
obvious if one identifies ±60° directions in Fig.[4]with y- and 
z-axes). Reducing the cubic lattice model to its F-model ver- 
sion by requiring that all vertex weights are invariant under 
Oh (the point symmetry group of the cubic lattice), we are left 
with only two vertex weights a and c (in notations of Fig. |4] 
b = a) and three switch weights a, p and £ (y = a and 8 = p). 
It then follows from Eqs. @ that 2a = c + 4p + 2£. If 2a > c, 
the model is in the DO phase. At 2a = c we have p = £ = 
meaning that loops must turn at every step, this is a 3D analog 
of non-self-intersecting loops. Every loop is perfectly AFE- 
ordered, and as we have seen, on a 2D square lattice this sig- 
nifies a transition to the AFE phase. There is also a general 
statement about self-intersections being a relevant perturba- 
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tion in a similar type of 2D loop models l37ll suggesting that 
the non-intersecting fully-packed loops characterise an unsta- 
ble fixed point. So once again, it is suggestive that 2a = c 
is a transition point. In 3D, however, the region with 2a < c 
is analogous to Region III of a triangular lattice model. In 
particular, there are zero modes in some of its ground states 
similar to that shown in in Fig.|5J; (naturally, these modes cor- 
responds to the loops of even length now). Thus, sorting out 
the properties of Region III of the F-model on the triangular 
lattice might shed some light on the properties of the cubic 
F-model. 



VI. CONCLUSION 

To conclude, we introduced a class of intersecting polymer 
loop models with the loop fugacity n. For n = 2, such a model 
is equivalent to an ice-type vertex model in its DO critical 
phase. We conjecture that the boundaries of the DO phase are 
encountered at the parameter values which force one or more 
switch weights to vanish. The onset of a long-range order in a 
vertex model is often, but not always equivalent to the forma- 
tion of cross-links between polymer loops that results in a vul- 
canisation transition. This correspondence is not strictly 2D, 
although in this paper we concentrated mostly on 2D exam- 
ples to benefit from the additional information available from 
the exact solutions. We note that all known transition points in 
the ice-rule models satisfy our criterion. As an additional ben- 
efit, we can understand transitions on irregular lattices [52]. 

However, despite its suggestive nature, it remains a big 
challenge to turn our criterion into a rigorous statement. 
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APPENDIX A: UNIFORM POLARISATION AND 
SUSCEPTIBILITY. 

Consider an L x L cylinder (or its generalised d- 
dimensional version - an L x L x . . . x L lattice with periodic 
boundary conditions imposed in d — 1 dimensions). The un- 
derlying lattice needs not be square or even regular as long 
as appropriate periodic boundary conditions can be imposed. 
The boundary condition in the remaining direction, i.e. along 
the cylinder axis is left free. We now apply a uniform electric 
field E along the axis. The coupling to this field is assumed to 
be of the standard form: AH = — E • dy where dy is a vec- 
tor in the direction of the arrow along the (ij) bond with the 
magnitude proportional to the length of the bond (for a regular 
lattice it is sufficient to choose all |dy | = 1). Because of the 
last condition, £/,-a 6 q E • dy =0 along any closed loop on the 



the vertex model in the form similar to Eqs. (^|2j: 
Z = E fl Ctf^W (2coshp£L) m ^ 

S k=\ 



(Al) 



where £ (s ) is the number of closed loops while m (s ) is the 
number of paths connecting the top and the bottom boundaries 
- see Fig.|6] The switch weights in Eq. JA1I correspond to the 
vertex weights in the absence of an external field. 




E 



FIG. 6: Only paths connecting the top and the bottom of the cylin- 
der can be polarised by the electric field parallel to its axis. The 
contribution of each such paths is proportional to the cylinder height 
regardless of the path's geometry. 



The uniform polarisation per site is then 



1 3lnZ 
N$ dE 



L l - d (m) tanhp£L >Z/~ d p£(m) (A2) 



2-dr 



for small field. Here (m) is the expected number of top-to- 
bottom paths. The uniform susceptibility is then 



dp 
dE 



,2-d 



m). 



(A3) 



£=0 



An extensive study of the square lattice six-vertex model in 
a uniform electric field presented in Ref. 1 15] suggests that this 
model has a finite non-zero susceptibility to such field every- 
where in the DO region. According to Eq. JA3t this implies 
that (m) remains finite in the thermodynamic limit. The fact 
that a probability of a path connecting the opposite sides of 
anLxL sample remains finite as its size grows is consistent 
with the r~ 2 decay of the loop correlation function. This is 
indeed the case for the points soluble via a free fermion map- 
ping 1 36, 38]. On the other hand, a mapping of the six-vertex 
model to an effective Gaussian model (via a restricted solid- 
on-solid (RSOS) model) and a Coulomb gas 13911 leads one 
to think that the exponent characterising the power law de- 
cay of correlations continuously changes throughout the DO 
phase. We notice, however, that the intersecting loops defined 
by the switch model are not the same as the non-intersecting 
loops running along the steps of equal height in the RSOS 
model which are described by the "magnetic" exponents in the 
Coulomb gas language. It must be the case then that the proper 
correlation function for our loops is described by a combina- 
tion of magnetic and electric operators "conspiring" to keep 



lattice. This allows one to evaluate the partition function of its r long-distance behavior throughout the DO region. 
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We finally remark that it has been conjectured in Ref. |40] 
that an arrow-arrow correlation function for a generic ice-rule 
model decays as r~ d (although, it seems, the basis for this 
claim is an exact result of 1 38] whose scope is limited to the 



free-fermion line of the 2D six-vertex model). This conjec- 
ture appears natural in light of Eq. jA3t : such behaviour is 
necessary for (m) ~ L d ~ 2 which in turn is required for a finite 
non-zero susceptibility to a uniform electric field. 
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